Scalar perturbation of the viscosity dark fluid cosmological model 
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Abstract 

A general equation of state is used to model unified dark matter and dark energy (dark fluid), and it has been 
proved that this model is equivalent to a single fluid with time-dependent bulk viscosity. In this paper, we investigate 
scalar perturbation of this viscosity dark fluid model. For particular parameter selection, we find that perturbation 
quantity can be obtained exactly in the future universe. We numerically solve the perturbation evolution equations, 
and compare the results with those of ACDM model. Gravitational potential and the density perturbation of the 
model studied here have the similar behavior with the standard model, though there exists significant value differences 
in the late universe. 
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I. INTRODUCTION 



Astrophysics and cosmology observations in recent years delineate the cosmological picture on its con- 
stituents more and more accurately, i.e, the precision cosmology era comes. Except for the long standing 
puzzling dark matter component, an unknown cosmic matter-energy constituent refereed to the so called dark 
energy may also exist that accelerates our universe expansion now, which contradicts with our traditional 
attitude on the behavior of conventional matter but it is eventually confirmed by recent observations like 
SNe la [ll 0] and CMB observations The cosmological dark sector, often divided as the mysterious DM 
and DE sectors respectively, takes around 95% of total energy budget of our universe. The concord ACDM 
model could be consistent with most of global astrophysics observational results. But the introduction of 
the cosmological constant simultaneously results in the yet to answer problems related directly to how to 
understand the fundamental physics theory, like the fine-tuning and the coincidence problems respectively. 
At the same time, "most" of course does not equal to "all" , some astrophysics problems still need be clarified 
and solved in the framework of the ACDM model Q , such as the the core singularities of the cold matter halo 
profiles. With the aim to understand the cosmic acceleration or dark energy phenomena, many theoretical 
models have been proposed, like the scalar field models and the modified Einstein gravity models [H| H 0| • 

Due to the limited scope of our experimental and observational tools, we have not yet been able to under- 
stand the "dark" nature and to detect the origin of DM and DE. Purely gravitational probes can not provide 
enough information to differentiate these two kinds of mysterious constitutions either. Therefore, from the 
phenomenological and practical point of view, a single (unified DM with DE) fluid description may be more 
plausible at least in the cosmic evolution description, which utilizes a single equation of state to model the 
dark matter and dark energy contributions together [13|. Generally, such models has a 



non-constant equation of state (EoS), which reflects both its dynamical and thermodynamics characters. The 
density dependent equation of state is widely investigated, such as the famous Chaplygin gas and generalized 
Chaplygin gas models, which assume an EoS form like the p = —A/p°' where the a is a model parameter. 
Another practical method to modify the EoS is by the introduction of cosmic viscosity media contribution 
which replaces the simplest perfect fluid EoS on a more physical and realistic basis. In the homogeneous 
and isotropic Friedmann-Robertson- Walker frame, only a bulk viscosity term which behaves as an additive 
pressure contribution can mimic both the two dark components and their coupling effects by playing a main 
role to influence the cosmic evolution. Different forms of the viscosity coefficients have been proposed like the 



density p dependent 14| or the redshift z dependent 
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In this letter, we will proceed our effort on the investigation of the viscosity dark fluid model to study the 
scalar perturbation of this model, for perturbation analysis can provide us a powerful tool to differentiate and 
constrain cosmology models finely as the calculation of perturbation quantities links the theoretical models 
with more plentiful and precise observations, like the cosmic microwave background (CMB) and large scale 
structure observation. There have been some researches on the perturbation evolution of the viscosity models 
[3 [l3 , by which we know that after corresponding model parameters chosen properly, the Chaplygin gas 
formulation can be viewed as a special case of the density dependent viscosity model. In the non-perturbative 
(zero order) level, the Chaplygin gas model can be exactly solved and fit the observational data well. But it 
has been found that in the perturbation level, there exist some unacceptable behaviors, like the blow up of 
density perturbation evolution and other peculiar behaviors [0 [l^ - One motivation to build other kinds of 
the viscosity models is to overcome these difficulties the Chaplygin gas models possess. Here, we will consider 
a time-dependent viscosity coefficient model, which is equivalent to the introduction of a general Equation of 
state(EoS) 0. The general EoS is 

p = (7 - l)p + Po + whH + wh2H^ + WdHH. (1) 
In the background, this model can fit the current astrophysics observational datasets consistently. We derive 
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its perturbation equations that govern the evolution of gravitational potential and density perturbation below. 
We numerically solve the perturbation equation, and compare it with that of conventional ACDM model and 
the Chaplygin gas model finding that the dark fluid model behaves well in different scales. Though there exists 
some value difference between the ACDM model and the dark fluid model in the late time evolution, their 
gravitational potential and density contrast shape and evolution behavior are similar by plotting respectively. 

This paper is organized as follows: In Sec. II, we summarize the calculations of scalar perturbation, and give 
the general evolution equation of the gravitational potential. In Sec. Ill, we briefly review the background 
evolution of the dark fluid model. In Sec. IV, we discuss the perturbation evolution of the dark fluid model. 
In Sec. V, we numerically solve the perturbation equation and compare it with other models. Finally, we 
present the conclusions in the last section. 



II. CALCULATIONS OF SCALAR PERTURBATION 



In this paper, we choose Newtonian gauge to calculate the scalar perturbation 

ds^ = -(1 + 2(f))df + a{tfSij{l ~ 2ij)dx'dx^. (2) 

If making the assumption here that there is no contribution from anisotropy inertia, it concludes that (p — ip. 

Generally, Einstein field equation with perturbed metric takes the form(for simplicity, we set k = \ here- 
after.) [aH [221 

d • d^ , 1 1 , 

0+^0 = -^(/O+p)<5m, (4) 

?i + 3-0+(2- + i^)0 = i(Jp + 3Jp), (5) 
a a 6 a'^ D 

where 5p and 5p are first order perturbation to zero-order cosmic density p and pressure p respectively. 
Perturbation to velocity of cosmic fiuid Sui is decomposed as 5ui = V^^u -f 5u'^. 5u is the scalar velocity 
potential. 5u[ is a divergenceless vector, which we also assume here contributes no effect. V denotes gradient 
with respect to comoving coordinate. From Eq. ([3]) and we obtain a constrain on first-order perturbation 
quantity 4>, Sp and Su 

d 

2—(l)-dp + 3-ip + p)Su^0. (6) 
a-^ a 

Also momentum and energy conservation equation to first order in perturbation could be derived 

6p + dt [{p + p)Su] + -{p+ p)5u + (p + p)0 = 0, (7) 
a 

5p+—{5p + 5p)+^^[\{p + p)5u]-2,{p + p)^ = Q. (8) 
A generalized parameterized equation of state may have the form as 

p=(7-l)p + /(p;p;aO + ... (9) 



or 



p={^-l)p + g{H;H;Pi) + ..., (10) 
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where ai and /3i are model parameters. When the models are discussed in flat universe k = 0, two parame- 
terized function f{p) and g{H) are categorized in one class. The equation of state reduces to the perfect fluid 
case with EoS as p = (7 — l)p when the model parameters ai vanish. 

For a barotropic equation of state, adiabatic sound speed is defined as = hence the ratio of pressure 
and density perturbation is c^. Therefore we could eliminate 6p from Eq. ([3]) and (l5|) ,then obtain the equation 
governs the gravitational potential 

+ (4 + 3cl)U + [2^ + (1 + ScDi^r - c'X]^ = (11) 
If write metric perturbation as Fourier integral 

0(x,O- j d3qe''i-0,(t), (12) 

we can make the substitution — > , where q is the wave number. 



III. BACKGROUND EVOLUTION 

In this section, we briefly review the background evolution behavior of the viscosity dark fluid model. A 
general form of EoS investigated in [20] is 

p = {'^ - l)p + Po + whH + wh2H^ + WdHH (13) 

One can prove that this generally parameterized EoS can be effectively equivalent to a single fluid with a 
time-dependent bulk viscosity: 

C = Co + Ci- + C2-, (14) 

a a 

and three parameters in the viscosity coefficient correspond to EoS parameters as 

WH = -3Co, (15a) 
WH, = -3(Ci+C2), (15b) 
WdH = -3(2- (15c) 

Flat Friedmann-Robertson- Walker metric reads 

ds^ = -dt^ + a(i)2^yda;Mx^ (16) 

The energy-momentum tensor with modified EoS could be written as 

T^.^ pU^U,+pHf,,, (17) 

where p represents modified pressure, in model concerned it is Eq.([T]), and in comoving coordinate = 
(1, 0, 0, 0) . Due to the correspondence between this modified EoS model and viscosity model, pressure could 
also he p = p ~ (9, where 9 — U^^ — 3a/ a. 

Using the EoS above and Friedmann equation, the equation of scale factor a(t) evolution could be obtained 

a _ -(37-2)/2- (KV2)wH2 + (K^)wrfg .Q^2 ^ -{k'^)wh a ^ ~{k^/2)po 
a l + {K^)wdH ^a' l + (K^/2)wdHa 1 + {k"^ /2)wdH' 

After redefining model parameters, there will be a compact form of evolution equation, at the same time, this 
form is comparable to perfect fluid case 
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where 



7+ /3)wg2 , . 

7 = -I , . 2/o^ ' ^^^) 

1 _ -(^V2)u;h) 
Ti 1 + (KV2)w;<iH ' 

^ ^ -(AtV2)po 
Tl 1 + {K^/2)wdH ' 

T2 r| ■ 



(21) 



(22) 



(23) 



The solution of scale factor is 



«W = «o{^(l + 7^or-^)exp[^(i + ^)] + 

^(l-7^oT+|)exp[_i^(l-l)]}^/^r (24) 

From Friedmann equation, cosmic density evolution reads 

-(1 + 7^?oT - ^)(i + ^)exp(^^) - (1 - 70or + ^)(^ - ^ 



P(i) 



3k^7^ 



(1 + 7^07 - ^)exp(i^) + (1 - 70or + 



(25) 



The case above is for 7 7^ 0, when take the limit of 7, solution could be obtained 



a{t) = aoexp[(i0ori + J)(exp(^) - l) - ^i^^] , (26) 



and cosmic density evolution 

p^t) = -|[i0oexp(^) + ||(exp(^) - 1)]. (27) 

Using Friedmann equation, Eq. ([1]) could be converted into a form, r.h.s of which is only the function of 
density p 

P=(7"l)p-|-A/f-4, (28) 



where parameters are defined the same as l|20 )) - (|23|) . 

IV. EVOLUTION OF SCALAR PERTURBATION 
A. Evolution equation 

In adiabatic perturbation case, ratio of adiabatic density and pressure perturbation equals to the adiabatic 
sound speed 

^ = cl (29) 

OPa 
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In this paper we pay our attention on single fluid model and investigate perturbation in adiabatic region, 
that is, pressure perturbation is proportional to the density perturbation. For a barotropic EoS model, the 
proportional efficient between pressure and density perturbation is merely a function of density p. On the 
other hand, in the condition of single (dark) fluid, we assume that dark sector interacts with baryon matter 
and the least dominated radiation negligibly. 

With the barotropic EoS (|28)) . adiabatic sound speed is = 7 — 1 — From Eq. (ITT|) . we obtain 

the correspondent equation of single (dark) fluid model 

4>g + fk + 50, = 0, (30) 

where two coefficients are defined by 

f[H) = {Z^ + l)H-^ (31) 

and 



-i^-(7-l-^)^. 



giH:q) = 2H + 3jH'--H-{^-l-——)^. (32) 



Conveniently, one can decompose metric perturbation as (j)q{t) = v{t)p{t), therefore obtains a differential 
equation about vit) and pit) 

v + (2- + J)v + { )v = 0. (33) 

P P 



If we choose the function p{t) properly as 



1 



p(i)=exp(-- / dt'fit')), (34) 



then the damping term can be eliminated, so Eq. (|33| reduces to 

v-ilf+\f-9)v = 0. (35) 

which takes a harmonic oscillator form. For the short wave limit, we have v + g{q] H)v — and g ~ 
— (7 — 1 — 37^)^2- If WKB condition is fulfilled, we have the approximate solution 

« ~ g(g.ff)i/4 {^+^^P(' / ^i^9{q;H)) + c_exp( - i j dt^g{q;H))} (36) 



B. The future solution of the gravitational potential 



Eq. (|30|) is too complicated to solve exactly, so here we will consider a simpler asymptotic case and try to 
extract the solution in this limit. If we strict that parameter Ti should be negative to confirm that sound 
speed is real, then we could see from Eq. (19) in j2Q,] that cosmic density approaches a constant value 



3T] 

-2 



as t — > 00, and negative Ti also confirms the positive of energy. It concludes that the universe will enter a de 
Sitter period then: H — > Ha, where Ha is a constant. Also the scale factor evolves exponentially: 



a(0=e^-* Ha = J-^, (38) 
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Eq. (I55t becomes 

« - [J/a - 5a('?)]« = 0, (39) 

where /a = (87 + 1)H\ — ^ is the late time asymptotic constant of function f{H) and gA{q) = 'ilH\ — 
1 1 ^ 

jr-7?A — (7 — 1 — 3T\lTX')iii2 ^'^ function of g{H\ q) when Hubble parameter approaches constant. Then we 
can solve the differential equation and get 



v(t) - - \) + C2 J-,„(^^-^^)r(i + -^), (40) 

where two parameter m and n are defined for simplicity by 

" = \SI - ^Vl + Y^H^ = 1 [(37 - l)»Hj - ^1', (41) 

When the physical wavelength is much longer than the Hubble radius 

Xphys » ^ « H. (43) 
a{t) 

The solution in this long wave limit could be obtained directly from Eq. (|30p . Assuming the solution takes 
the form as 

0g(i)=e''*. (44) 

After inserting it into Eq. (pO|) . we get 

This quadratic equation has two dependent solution, therefore the solution of Eq. (PP)) reads 

— )t + caexp (-67iJA + 

This unstable solution exponentially increases or decreases in the infinite future, which dependents on the 
parameters. 

For the opposite limit, we consider large wave number solution with q 3> aiJ, and fix at some conformal 
wave number k — q/a. Then we have the simple solution in this limit 

V oc e^*, (47) 



52 + [(37 + 1)Ha - —]b + 3jHl - —Ha 0. (45) 



; = ciexp[(-//A + — r^)t] + c2exp[(-67i?A + —7Tr^)t] (46) 



where / = 7 — 1 — ^TiH ■ Together with the definition p4l) , we see this is an exponential decay solution for 
the gravitational potential, which is also unstable. 



V. NUMERICAL RESULTS 
A. Comparison models 

In the numerical results presented below, we also calculate perturbation in ACDM model and the Chaplygin 
gas model numerically as a comparison, so we present a short review here. 

(a) ACDM model: The cosniological constant does not contribute perturbation in the total energy density 



SpA = 0. (48) 
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The density perturbation comes from matter density, Sp — Spm- The cosmological constant plays a role in in- 
fluencing the background evolution of the universe, especially by Hubble parameter, so it will be imprinted into 
the evolution of matter perturbation. We could obtain the scale-independent gravity perturbation equation 
from Eqs. ^ and 

4>mq + 4-ff</)™g + {2H + ?>H^)(t>mq = 0- (49) 

In the matter dominated era 

= ^Prn = ^mOd-^ = ^t'^ ■ (50) 

Also the solution ipmq = (j^mqa + ^t^^^/^ is well-known, 
(b) Chaplygin gas model: It has the EoS form 

P = -- (51) 
P 

If we use this EoS to model single fluid universe, then get Hubble parameter 

H(z) = p{l + zf + (52) 



and the adiabatic sound speed 



Hence H and in Eq. are specified. 



A 



< = - (53) 
P 



B. Effective dark fluid 

1. Background 

For simplicity, we set parameter T2 — 00, which we call this case effective dark fluid. This means if p — 0, 
pressure p vanishes too, and there is not a cosmological constant like pressure contribution, which could be 
seen from EoS (pS)) . Hence Eq. becomes 

'--^'^i^f^—. (54) 
a I a ±1 a 



This differential equation about scale factor can be seen as a special case given by a general EoS [2J| [25| [26| 

p = -p- Ap°' - BH^f', (55) 

which can give 

-= ^ - +A- +M- +'^- 56 

a 2 a' ^a' 

Eq. (|54p can be converted as a differential equation of H(a) 
Its solution is 

H{a)/Ho = na-^"*/^ + {1 - n), (58) 
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where SI = 1 — ^yKHo' '"^"^^ have already set qq = 1. We use this Hubble parameter with 7 = 0.9 and 
7 = 1.2 to calculate the distance 



DLiz)^Hail 



H{z') 



(59) 



and the distance modulus. We compare it with the supernova data 23|, which is plotted in FIG. 1. We see 
that model with parameter 7 = 0.9 (blue) and 7 — 1.2 nearly can not be discriminated in the late time (small 
redshift, more data in this region has been obtained). There are some differences for larger redshift. But in 
the whole, both two parameters consist with the data well. 




FIG. 1: The relation between distance modulus and redshift. Black and Blue lines correspond to the theoretical 
calculation curves with 7 = 1.2 and 0.9 respectively. Both of them can fit data in an acceptable level. 



2. The gravitational potential and the density perturbation 

After inserting the Hubble parameter into Eq. (pO)) . we numerically solve this scale-dependent model. 
Results are illustrated in FIG. 2. As a comparison, ACDM and Chaplygin gas model are solved and plotted 
(green and red line respective) too. In the early time, the difference between ACDM and dark fluid model 
is tiny. Both of them behaves nearly as a constant. In the late time, two models give the same shape of the 
gravitational potential. Though these two models predicts the decay of potential, there exists a value contrast 
around 5% — 10%. The quantity of contrast dependents on the value of model parameter 7. Here we point 
that the dark fluid model predicts more similar potential as ACDM than the Chaplygin gas model. On the 
other hand, it can be qualitatively seen the level of scale-dependent of the dark fluid model. We plot different 
results with q — 0.005, 0.5 and 1.5. 

For the small scale, the perturbation evolution of the dark fluid model with 7 = 0.9 is significantly different. 
The gravitational potential decays much earlier, which contradicts with ACDM. Numerical results for 7 < 1 
indicate that this condition strongly influences the early evolution of perturbation quantity(we only plot 
7 = 0.9 case here.), which puts strict constrain on parameter region. Also, if 7 is bigger than 1, the positive 
of adiabatic sound speed could be easily fulfilled. 

From Eq. ([3]), density perturbation could be expressed as 

5p^~&H(})~&H'^(f) + 2^<i). (60) 

It tells us that once we have the solution of the gravitational potential we could get the information of the 
density contrast. During the early time (0 ~ 0) and in the very large scale {q <C ai?), the density contrast 
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FIG. 2: The gravitational potential for different wave number(g — 0.005, 0.5, 1.5) evolves with time. Black, blue, green 
and red lines represent the dark fluid model with 7 = 1.2 and 0.9, ACDM and the Chaplygin gas model respectively. 



and the gravitational potential is linked by the Hubble parameter Sp ~ —6H^(j). Always define 

where Friedmann equation is used in the last step. Numerical results of the density perturbation S for different 




Log,„a 



FIG. 3: Density perturbation with q = 0.005, 0.5 and 1.5 evolves with time. Black, blue, green and red lines represent 
the dark fluid model with 7 = 1.2 and 0.9, ACDM and the Chaplygin gas model respectively. 



q is plotted in FIG. 3. The numerical curves have the similar shape, but for modified models, the density 
perturbation is suppressed in the late time. For the large scale, the density perturbation evolution in the 
dark fluid model increases linearly in the late universe. 5 deviates from that of ACDM in the late time, and 
suppressed today for different scale. The values enhanced or depressed dependent on parameter 7 and scale. 



VI. CONCLUSION AND DISCUSSION 



In this paper, we investigate extensively the dark fluid model proposed in [20|, equivalently this model can 
be viewed as a single fluid with time-dependent bulk viscosity. Scale factor and density evolution can be 
exactly solved in this model. In the background, the dark fluid model can fit the supernova data acceptable. 
Our main task in this paper is to analysis the behavior of this model in the perturbation level. We derive 
equations govern the perturbation quantities. For the condition that Ti is smaller than 0, the universe will 
enter de Sitter phase in the t ^ 00 future. We solve exactly the gravitational equation in this condition and 
obtain the solution for both long and short wave case. Generally, the perturbation evolution equations are 
solved numerically. When compare the results with those of ACDM model, we find that 
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• In the early time and the large scale, both the gravitational potentials of two models behave as a 
constant. 

• Though the gravitational potentials of two models have similar behavior and shape, as can be seen form 
FIG. 3, there exists about 5% — 10% significant value difference in the late time. 

Perturbation analysis also provides constraint on model parameter. For different selection of parameter 7, 
both 7 < and 7 > can give consistent prediction curve of distance modulus, but numerical results indicate 
in 7 < case, the gravitational potential deviates from ACDM significantly from the early time. This result 
strongly constrain the selection region of 7. We suggest 7 > 0, which can also produce positive sound speed 
naturally. 
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